This is a first-pass exploration of some of the relationships among QTS data for the “Spiritual curiosity” interview.
Requested by Josh
prayfreqminHow many minutes each day do you do that on average throughout the week? [use real world examples, ask alone? With others? Probe away– Answer with interviewer judgement for answers]
godvoxaloudSome people say that they have heard God* speak out loud to them. Has this ever happened to you?
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godvoxaloud_num ~ prayfreqmin + country, family = "binomial",
data = d1)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.703 -1.027 -0.455 1.032 2.154
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.47926 0.13784 -3.477 0.000507 ***
prayfreqmin.L 0.73747 0.24453 3.016 0.002562 **
prayfreqmin.Q 0.54372 0.27178 2.001 0.045436 *
prayfreqmin.C -0.35794 0.29205 -1.226 0.220347
countryGH 0.97536 0.21961 4.441 8.94e-06 ***
countryTH 0.25817 0.26628 0.970 0.332274
countryCH 0.30449 0.24600 1.238 0.215809
countryVT 0.05578 0.25920 0.215 0.829609
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 446.17 on 328 degrees of freedom
Residual deviance: 385.93 on 321 degrees of freedom
AIC: 401.93
Number of Fisher Scoring iterations: 4
godviavisions_Some people say that they have had a vision from God*—they have a picture, but it is like they see it with their eyes. Has anything like that happened to you?_
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godviavisions_num ~ prayfreqmin + country, family = "binomial",
data = d2)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.8640 -1.0136 -0.5422 1.1679 2.0691
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.6210 0.1361 -4.562 5.06e-06 ***
prayfreqmin.L 0.9386 0.2470 3.800 0.000145 ***
prayfreqmin.Q 0.8468 0.2739 3.092 0.001989 **
prayfreqmin.C 0.3692 0.2926 1.262 0.207031
countryGH 1.0290 0.2239 4.597 4.29e-06 ***
countryTH 0.6839 0.2719 2.515 0.011894 *
countryCH -0.9077 0.2933 -3.095 0.001968 **
countryVT 0.1282 0.2570 0.499 0.617930
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 443.63 on 334 degrees of freedom
Residual deviance: 386.58 on 327 degrees of freedom
AIC: 402.58
Number of Fisher Scoring iterations: 4
godviabodyexperiencesSome people have particular experiences in your body that they associate with God* or spirit. Does that happen for you? [examples: warm hands, goosebumps, fire in the belly]
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godviabodyexperiences_num ~ prayfreqmin + country,
family = "binomial", data = d3)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.0197 -1.1520 0.6754 0.8918 1.3475
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 0.7061 0.1332 5.302 1.14e-07 ***
prayfreqmin.L 0.5249 0.2378 2.207 0.02731 *
prayfreqmin.Q 0.5157 0.2622 1.966 0.04926 *
prayfreqmin.C 0.2188 0.2779 0.788 0.43096
countryGH 0.5353 0.2375 2.254 0.02418 *
countryTH 0.8843 0.3073 2.877 0.00401 **
countryCH -0.4790 0.2488 -1.926 0.05415 .
countryVT -0.1303 0.2500 -0.521 0.60225
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 432.22 on 334 degrees of freedom
Residual deviance: 402.36 on 327 degrees of freedom
AIC: 418.36
Number of Fisher Scoring iterations: 4
Linear mixed model fit by REML ['lmerMod']
Formula: response ~ prayfreqmin + country + (1 | subject_name) + (1 |
question)
Data: d_agg1
REML criterion at convergence: 9516
Scaled residuals:
Min 1Q Median 3Q Max
-2.6638 -0.7863 -0.1873 0.8822 2.5579
Random effects:
Groups Name Variance Std.Dev.
subject_name (Intercept) 0.02090 0.1446
question (Intercept) 0.03517 0.1875
Residual 0.17972 0.4239
Number of obs: 7961, groups: subject_name, 344; question, 24
Fixed effects:
Estimate Std. Error t value
(Intercept) 0.459403 0.039550 11.616
prayfreqmin.L 0.117263 0.017788 6.592
prayfreqmin.Q 0.033800 0.019728 1.713
prayfreqmin.C 0.048443 0.021485 2.255
countryGH 0.024153 0.016975 1.423
countryTH 0.001581 0.020844 0.076
countryCH -0.086277 0.019517 -4.421
countryVT 0.155426 0.019417 8.005
Correlation of Fixed Effects:
(Intr) pryf.L pryf.Q pryf.C cntrGH cntrTH cntrCH
prayfrqmn.L 0.035
prayfrqmn.Q -0.066 0.174
prayfrqmn.C 0.024 -0.164 -0.026
countryGH -0.045 -0.099 0.114 0.118
countryTH 0.034 0.240 0.040 -0.020 -0.259
countryCH 0.014 -0.076 -0.095 -0.034 -0.227 -0.305
countryVT 0.004 -0.045 0.058 -0.022 -0.209 -0.299 -0.278
Linear mixed model fit by REML ['lmerMod']
Formula: response ~ prayfreqmin * country + (1 | subject_name) + (1 |
question)
Data: d_agg1
REML criterion at convergence: 9563.5
Scaled residuals:
Min 1Q Median 3Q Max
-2.6672 -0.7887 -0.1796 0.8787 2.5708
Random effects:
Groups Name Variance Std.Dev.
subject_name (Intercept) 0.02118 0.1455
question (Intercept) 0.03517 0.1875
Residual 0.17973 0.4239
Number of obs: 7961, groups: subject_name, 344; question, 24
Fixed effects:
Estimate Std. Error t value
(Intercept) 0.458378 0.039986 11.463
prayfreqmin.L 0.108947 0.022863 4.765
prayfreqmin.Q 0.026234 0.023102 1.136
prayfreqmin.C 0.042343 0.023338 1.814
countryGH 0.027652 0.018607 1.486
countryTH -0.012737 0.030968 -0.411
countryCH -0.075852 0.022478 -3.374
countryVT 0.153817 0.020838 7.381
prayfreqmin.L:countryGH 0.006461 0.035854 0.180
prayfreqmin.Q:countryGH 0.003627 0.037213 0.097
prayfreqmin.C:countryGH 0.023226 0.038524 0.603
prayfreqmin.L:countryTH -0.039139 0.069967 -0.559
prayfreqmin.Q:countryTH -0.016074 0.061933 -0.260
prayfreqmin.C:countryTH 0.009706 0.052689 0.184
prayfreqmin.L:countryCH -0.029814 0.038224 -0.780
prayfreqmin.Q:countryCH -0.039246 0.044956 -0.873
prayfreqmin.C:countryCH 0.018775 0.050805 0.370
prayfreqmin.L:countryVT -0.012807 0.039717 -0.322
prayfreqmin.Q:countryVT 0.036971 0.041676 0.887
prayfreqmin.C:countryVT 0.035715 0.043548 0.820
Correlation matrix not shown by default, as p = 20 > 12.
Use print(x, correlation=TRUE) or
vcov(x) if you need it
refitting model(s) with ML (instead of REML)
Data: d_agg1
Models:
r_agg1: response ~ prayfreqmin + country + (1 | subject_name) + (1 |
r_agg1: question)
r_agg2: response ~ prayfreqmin * country + (1 | subject_name) + (1 |
r_agg2: question)
Df AIC BIC logLik deviance Chisq Chi Df Pr(>Chisq)
r_agg1 11 9489.8 9566.6 -4733.9 9467.8
r_agg2 23 9504.7 9665.3 -4729.3 9458.7 9.134 12 0.6914
Ignoring unknown aesthetics: y
Requested by Emily, John
selfunsuregodrealHas there been a time when you yourself wondered whether God* was real?
godvoxaloudSome people say that they have heard God* speak out loud to them. Has this ever happened to you?
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godvoxaloud_num ~ selfunsuregodreal + country,
family = "binomial", data = d4)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.4270 -1.0266 -0.5155 0.9986 2.0415
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.47400 0.12983 -3.651 0.000261 ***
selfunsuregodrealY -0.06683 0.13151 -0.508 0.611315
countryGH 0.97704 0.21634 4.516 6.29e-06 ***
countryTH 0.01075 0.26818 0.040 0.968028
countryCH 0.38097 0.24572 1.550 0.121032
countryVT 0.04141 0.26124 0.159 0.874046
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 426.65 on 313 degrees of freedom
Residual deviance: 381.81 on 308 degrees of freedom
AIC: 393.81
Number of Fisher Scoring iterations: 4
godviavisions_Some people say that they have had a vision from God*—they have a picture, but it is like they see it with their eyes. Has anything like that happened to you?_
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godviavisions_num ~ selfunsuregodreal + country,
family = "binomial", data = d5)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.3530 -0.9346 -0.6884 1.0629 1.8214
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.61222 0.12652 -4.839 1.30e-06 ***
selfunsuregodrealY 0.06415 0.13016 0.493 0.6221
countryGH 0.95182 0.21279 4.473 7.71e-06 ***
countryTH 0.33395 0.26614 1.255 0.2096
countryCH -0.58890 0.27649 -2.130 0.0332 *
countryVT 0.07429 0.25358 0.293 0.7696
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 422.46 on 318 degrees of freedom
Residual deviance: 391.05 on 313 degrees of freedom
AIC: 403.05
Number of Fisher Scoring iterations: 4
godviabodyexperiencesSome people have particular experiences in your body that they associate with God* or spirit. Does that happen for you? [examples: warm hands, goosebumps, fire in the belly]
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godviabodyexperiences_num ~ selfunsuregodreal +
country, family = "binomial", data = d6)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.8606 -1.2321 0.6274 0.9309 1.2922
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 0.7056 0.1277 5.523 3.33e-08 ***
selfunsuregodrealY 0.2421 0.1317 1.839 0.0659 .
countryGH 0.5778 0.2363 2.446 0.0145 *
countryTH 0.5884 0.3089 1.905 0.0568 .
countryCH -0.3357 0.2470 -1.359 0.1742
countryVT -0.1013 0.2519 -0.402 0.6876
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 411.83 on 319 degrees of freedom
Residual deviance: 390.07 on 314 degrees of freedom
AIC: 402.07
Number of Fisher Scoring iterations: 4
morequesmoreanswrDo you think that the more spiritually mature you become, you will discover more questions or more answers?
godvoxaloudSome people say that they have heard God* speak out loud to them. Has this ever happened to you?
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godvoxaloud_num ~ morequesmoreanswr + country,
family = "binomial", data = d7)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.4044 -0.9984 -0.4288 0.9900 2.2055
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.65574 0.16797 -3.904 9.46e-05 ***
morequesmoreanswrQ -0.03057 0.15819 -0.193 0.847
countryGH 1.14447 0.23706 4.828 1.38e-06 ***
countryTH 0.18839 0.31250 0.603 0.547
countryCH 0.28946 0.33757 0.857 0.391
countryVT 0.03154 0.32558 0.097 0.923
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 325.90 on 241 degrees of freedom
Residual deviance: 280.75 on 236 degrees of freedom
AIC: 292.75
Number of Fisher Scoring iterations: 4
godviavisions_Some people say that they have had a vision from God*—they have a picture, but it is like they see it with their eyes. Has anything like that happened to you?_
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godviavisions_num ~ morequesmoreanswr + country,
family = "binomial", data = d8)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.388 -1.039 -0.616 1.127 2.026
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.6647 0.1687 -3.939 8.18e-05 ***
morequesmoreanswrQ 0.1811 0.1525 1.188 0.2349
countryGH 0.9661 0.2353 4.106 4.03e-05 ***
countryTH 0.5124 0.3171 1.616 0.1061
countryCH -1.0697 0.4511 -2.372 0.0177 *
countryVT 0.3111 0.3162 0.984 0.3251
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 325.22 on 242 degrees of freedom
Residual deviance: 295.17 on 237 degrees of freedom
AIC: 307.17
Number of Fisher Scoring iterations: 4
godviabodyexperiencesSome people have particular experiences in your body that they associate with God* or spirit. Does that happen for you? [examples: warm hands, goosebumps, fire in the belly]
Ignoring unknown aesthetics: fill
Ignoring unknown aesthetics: fill
Call:
glm(formula = godviabodyexperiences_num ~ morequesmoreanswr +
country, family = "binomial", data = d9)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.9152 -1.2161 0.7310 0.9641 1.1674
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 0.71079 0.15975 4.449 8.62e-06 ***
morequesmoreanswrQ 0.03348 0.15584 0.215 0.8299
countryGH 0.43886 0.24850 1.766 0.0774 .
countryTH 0.91566 0.38084 2.404 0.0162 *
countryCH -0.54853 0.33456 -1.640 0.1011
countryVT -0.15232 0.30627 -0.497 0.6190
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 310.96 on 244 degrees of freedom
Residual deviance: 294.81 on 239 degrees of freedom
AIC: 306.81
Number of Fisher Scoring iterations: 4
Here are other things on our to-do list: